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QUASI -ONE-DIMENSIONAL COMPRESSIBLE FLOW 
ACROSS FACE SEALS AND NARROW SLOTS 
I - ANALYSIS 

by John Zuk, Lawrence P. Ludwig, and Robert L. Johnson 

Lewis Research Center 

SUMMARY 

An analysis is presented for compressible fluid flow across gas film shaft face 
seals and narrow slots. The flow is analyzed using a quasi- one- dimensional integral 
method which includes fluid inertia, viscous friction, and entrance losses. Viscous 
friction is accounted for by a mean friction factor. The model is valid for both laminar 
and turbulent flow. Subsonic and choked flow conditions can be analyzed. 

A seal radial width to film thickness geometric parameter is used. For large val- 
ues of this parameter and for low Mach number, the results are identical to the solution 
for viscous subsonic flow; for small values, the results correspond to orifice flow equa- 
tions. Results agree with radial flow experiments. 

Results show that a parallel film can have a positive film stiffness under choked 
flow conditions. 


INTRODUCTION 

Shaft seal systems in advanced aircraft turbine engines will be operated at speeds, 
temperatures, and pressures higher than shaft seals currently used. Conventional face 
seals presently used in gas turbine engines are limited to sliding velocities of about 
110 meters per second (350 ft/sec), pressures of about 86 newtons per square centime- 
ter (125 lb/in. ^), and gas temperatures of 700 K (800° F) (ref. 1). Advanced engines, 
however, will require seals to operate to speeds of 150 meters per second (500 ft/sec) 
(ref. 2), pressures to 340 new tons per square centimeter (500 lb/in. ), and tempera- 
tures to 980 K (1300° F) (ref. 3). For face seals operating at these conditions, a posi- 
tive face separation (no rubbing contact) will be required to achieve long life and relia- 
bility. A promising approach is the use of gas film seals with self-acting lift pads (see 
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Figure 1. - Pressure-balanced face seal with self-acting lift pads (added for axial film stiffness). 

fig. 1, refs. 4 and 5, and Seal Description and Models). 

Since the seals must be pressure balanced, a proper balance of the opening and 
closing forces must be achieved with a leakage gap that has tolerable mass leakage. 

The gap must be small enough so that the leakage is minimal, but it must be large 
enough so that power dissipation, due to shear in the film, and face deformation are 
tolerable. Thus the design of the sealing gap is vital to seal performance and pressure 
distribution in the gap and mass leakage through the gap must be analyzed. In this re- 
port only the sealing dam portion of the seal is analyzed. 

The classical viscous, isothermal, subsonic, compressible flow analysis for this 
problem is well known (e.g. , see Gross, ref. 6). The pressure distribution and mass 
leakage have been calculated for the parallel film hydrostatic case. Carothers (ref. 7) 
has conducted compressible flow experiments in thin films of air which are flowing 
radially between parallel plates. The pressure distribution was found for both subsonic 
and supersonic axial entrance flows. Carothers’ geometry, however, was representa- 
tive of externally pressurized gas bearings rather than seals. Grinell (ref. 8) has the- 
oretically and experimentally investigated compressible flow in thin passages and has 
shown agreement between theory and experiment. Mueller (ref. 9) included compres- 
sibility for the case where there was a restricted bypass orifice flow into the parallel 
gap, but his analysis did not consider conditions which would yield choked flow in the 
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gap. The textbook on mechanical face seals by Mayer (ref. 10) virtually excluded face 
seal operation under conditions where fluid compressibility is important. 

The objective of this report is to present a mathematical analysis of compressible 
fluid flow across shaft face seals and narrow slots. The analysis includes fluid inertia, 
viscous friction, and entrance losses and thus describes both subsonic and choked flow 
conditions. For specified pressure ratios, film thicknesses, and friction factor varia- 
tions with Reynolds number, which are used as parametric input, the output variables 
include mass and volume leakage flow rates, force, center of pressure, and distribu- 
tions of pressure, density, velocity, temperature, and Mach number. 


Seal Description and Models 

A conventional face seal which is pressure (force) balanced is shown in figure 2. 

In this seal, the pressure drop occurs across a narrowly spaced sealing dam, and the 
force due to this pressure drop is balanced by a hydrostatic closing force and spring 
force. This configuration, however, has an inherent problem. For parallel sealing 
dam surfaces, the force caused by the pressure drop across the sealing dam is inde- 
pendent of film thickness for low Mach numbers; hence, there is no way of maintaining 
preselected film thickness which will allow tolerable leakage and still have noncontact 
operation. Since the force is independent of film thickness, the design also lacks axial 
film stiffness for sufficient dynamic tracking of the stationary nosepiece with the rotat- 
ing seal seat. In order to operate adequately, the seal nosepiece must follow the seal 
seat surface under all conditions without surface contact or excessive increase in film 
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thickness, which would yield high leakage. Some of these conditions are axial runout, 
misalinement, thermal distortion, coning, and dishing. 

A promising method of maintaining a preselected film thickness and achieving axial 
film stiffness is to add a self-acting gas bearing, such as shrouded Rayleigh step lift 
pads, to the conventional rubbing contact pressure- balanced face seal (see fig. 1 for a 
sketch and refs. 4 and 5 for details of this concept, design, and test results). Both the 
sealing dam force due to the pressure drop across the sealing dam and the self-acting 
lift pad (gas bearing) force are balanced by the hydrostatic and spring closing forces. 
The gas bearing has a desirable characteristic whereby the force increases with de- 
creasing film thickness. If the seal is perturbed in such a way as to decrease the gap, 
the additional force generated by the gas bearing will open the gap to the original equili- 
brium position. In a similar manner, if the gap becomes larger, the gas bearing force 
decreases, and the closing force will cause the seal gap to return to the equilibrium 
position. 


Some Models and Their Limitations 

For subsonic flow, variable cross-sectional area, and Mach number < l/^y (all 
symbols are defined in appendix A), the analysis of reference 11 can be used. This dif- 
ferential analysis finds a solution for the case where the viscous friction is balanced by 
the pressure drop and the fluid inertia is negligible. For this condition the entrance 
pressure drop is small; hence, Pq = P^ and Pg = Pg (see fig. 3). The flow is nearly 
isothermal for the film thickness range studied. The analysis of reference 11 yields the 
classical cubic dependence of mass flow on film thickness. Also, this analysis enables 
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Figure 3. - Model and notation of sealing faces including control volume for 
quasi-one-dimensional flow. 
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small tilts of the sealing dam surfaces to be studied. These small tilts simulate seal 
face deformations which can occur due to thermal, centrifugal, etc., effects. 

For flow approaching Mach 1 (i.e. , choked flow), the inertia force neglected in 
reference 11 becomes important. Also, the flow behaves more as an adiabatic flow than 
an isothermal flow. The equations (partial differential equations) are relatively complex 
because of the nonlinearity of the inertia terms; thus an approximate analytical model 
must be used. In the approximate model presented here, the flow is analyzed as a 
quasi- one- dimensional flow using a control volume integral analysis. 

Since in gas film face seal applications there is seal seat rotation relative to the 
seal nosepiece, a question arises as to the validity of neglecting the centripetal inertia 
force which would affect the radial flow. This question was examined in reference 13. 

It was found that the circumferential rotational flow can be uncoupled from the radial 
pressure flow when the ratio of rotational velocity to pressure flow velocity is less than 
l/\/Re(h/AR). Rotational effects are important for calculating the viscous shear and 
power dissipation. 

It sh~ -d oe nointed out that the flow in a gas film face seal is qualitatively quite 
different from the flow in an externally pressurized thrust bearing. Both operate with 
very small film thicknesses, usually less than 0. 0025 centimeter (1 mil). However, 
gas film face seals are characterized by a small (Rg - R^)/Rj ratio. Also, the inner 
cavity is a large supply reservoir. Hence, in gas film seals the reservoir pressure 
does not vary with film thickness. On the other hand, externally pressurized gas thrust 
bearings are characterized by a large (Rg - Rj)/R^ ratio. (The large surface area is 
necessary for high load capacity. ) In order to maintain positive film stiffness, compen- 
sating inlet flow restrictors are used. Because of these restrictors, the inlet pressure 
to the bearing varies significantly with film thickness. As a result of these differences, 
the flow regimes for the two cases are quite different. 


QUASI -ONE- DIMENSIONAL FLOW ANALYSIS 

As stated in the previous section, when the radial flow is close to choked or is 
choked at the exit, the viscous flow analysis in reference 11 is no longer valid. Fluid 
inertia, neglected in that analysis, becomes important. The rigorous inclusion of iner- 
tial terms necessitates the solution of nonlinear partial differential equations and is 
therefore not practical for engineering calculations; hence an approximate analytical 
model will be formulated which will be especially useful for calculating flow conditions 
near and at choking. 

The geometry is shown in figure 3. Two parallel, coaxial, circular rings are sep- 
arated by a very narrow gap. Rotational effects on the radial flow are negligible for 
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most applications (ref. 13). A pressure differential exists between the rings’ inner and 
outer radii. The cavities on either side (i.d. and o.d. ) of the sealing dam are assumed 
to be constant pressure reservoirs. For subsonic flow, reservoir conditions and the 
exit ambient pressure are specified. For choked flow, reservoir conditions and an exit 
Mach number of one are specified. 

The analysis can be separated into two parts, which can then be considered sep- 
arately. One part is an analysis of the seal passage itself, in which the flow is assumed 
to behave as a constant area, adiabatic flow with friction, while the other part is an 
analysis of the entrance flow. The analysis for these two regions is described first fol- 
lowed by a discussion of the iterative procedure for the complete solution. 


Constant Area Adiabatic Flow With Friction 

It is assumed that the flow in the seal leakage flow region behaves as a constant 
area adiabatic flow with friction. A quasi- one- dimensional approximation is made 
wherein it is assumed that the flow properties can be described in terms of their cross- 
sectional averages. 

The following assumptions have been made in the analysis: 

(1) The area expansion due to radius increase is neglected. (In most mainshaft face 
seals the i. d. /o. d. is about 0. 98. ) 

(2) The flow is adiabatic. 

(3 ) No shaft work is done on or by the system . 

(4) No potential energy differences are present such as caused by elevation differ- 
ences. 

(5) The fluid behaves as a perfect gas. 

(6) The sealing surfaces are parallel. 

With these assumptions, the flow is commonly known as Fanno line flow (ref. 15). This 
analysis is similar to that used by Grinnell (ref. 8) for flow in thin passages. The gov- 
erning equations when area changes are neglected are as follows: 

Conservation of mass : 


M - puA = constant 


which reduces to 


dp + 1 du 2 = Q 
P 2 u 2 


( 1 ) 


( 2 ) 
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Conservation of energy: 


or 


e 



P 2 


= constant 


i + — - = constant 
2 


where i is the specific enthalpy. This can be written as 


di = c dT = - - du 2 

P 2 


This equation can be reduced to 


dT, (y - 1) 
T 2 



where the Mach number M = u/-y/y (RT. 
Equation of state: 


P = p(RT 


or 


dP = dp + dT 
P p T 


Conservation of momentum: 


Now the hydraulic diameter and Fanning friction factor are introduced: 


D = 


4A 


dA 


w 


dX 


2 

It is assumed that the viscous effect can be represented by a mean friction factor. The 
mean friction factor depends on Reynolds number only and is a slowly varying function 
of Reynolds number for many flows, geometries, and conditions. 

The set of equations ((2), (5), (7), and (8)) can be combined into a single equation in 
terms of the Mach number alone. Details of obtaining this equation can be found in 
appendix B. The result is 


«dX=. 
D 


(1 - M z ) dM z 
yM 4 ^1 + M 2 ^ 


(9) 


This equation can now be integrated from some location X (where the Mach num- 

a 

ber is M a ) to the location X^ (where the Mach number is M^). By the use of partial 
fractions, the right side of equation (9) can be integrated to 


where 


»M, 


yM 


r M„ 


2 

i 1 " M } dM 2 - B(M ) - B(M b ) 


B(M)=ii^M%X±iin i- 


y±A m 2 


r m 2 2 ^ l+iiiii 2 


do) 


(ii) 
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Since the variation of f with X is usually unknown, the integral of the left side of 
equation (9) cannot be evaluated exactly. However, if f is replaced by an appropriate 
mean value f , the integration may be performed: 


_4 

D 



(x b - X a ) 


The integration of equation (9) yields 

f ( x b " X a> = B ( M a) ' B < M b> 


( 12 ) 


(13) 


In the computer calculation, it is convenient to define a length 

L N = X max “ X N 


such that 


4f 

D 


(X. 


max 


X N>=‘ 


4f L 


'N 


D 



D 


- B(Mjq-) 


(14) 

(15) 

(16) 


(Note that B(M = 1) = 0. ) The length L N is the length of the seal (measured from X N ) 
for which the condition at the exit is sonic (choked). It follows, if and Mg are the 
Mach numbers at the inlet and exit of the seal, respectively, that 

AR = Lj - Lg (17) 


See figure 4. 

The reservoir pressure and exit ambient static pressure are known. It is desirable 
to find the pressure at any point in the flow leakage passage as a function of Mach num- 
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Figure 4. - Lengths and stations used in the analysis. (Subsonic case 
is shown. ) 


ber. The flow is adiabatic. Thus the stagnation temperature is constant everywhere in 
the seal; that is, 


T 


01 


- T 


02 


(18) 


or (ref. 15, p. 80) 


T i (» + ^ M ?) = t 2 (* + *7* M i) <19) 

From this relation, the continuity equation, the equation of state, and the definition of 
the Mach numbers, one obtains 



( 20 ) 


Entrance Flow 


If it is assumed that the flow in the entrance region is isentropic and adiabatic, the 
analysis is straightforward. The adiabatic energy equation is 


i 0 



( 21 ) 
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which can be written as 



(y ~ 1)M^ 

2 


( 22 ) 


By use of the isentropic relation, P/p^ = constant, and the equation of state, other inlet 
variables can be found knowing the sealed high pressure reservoir values; for example, 



It is well known that entrance flows do not behave as isentropic flows but that an 
additional pressure drop is present due to viscous friction, turning losses, etc. A prac- 
tical way of accounting for the entrance loss is to introduce an empirically determined 
velocity loss coefficient C-^ (see appendix C). Hence, equation (22) becomes 



(r- DMf 
1 + 1 



while equation (23) becomes 


P l = 


1 + 


(y - DM| 


2C( 


y/(y-i) 


(25) 


(26) 


Under the entrance conditions discussed in appendix C, the entrance velocity loss 
coefficient C L is related to the head loss coefficient k' commonly used in hydraulics 
by the relation 
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- 



(27) 


k* = 



1 


Iteration Procedure 


The procedure for solving the previous equations is an iterative one. The solution 
is found using a computer program which appears in reference 12. In the iterative pro- 
cedure, the first step is to solve for the minimum film thickness for which the flow is 
choked. If the actual film thickness is less than this minimum value, the speed of the 
flow is less than sonic at the exit and the exit condition is Pg = Pg. If the actual film 
thickness is greater than the minimum for choked flow, the exit flow is sonic and the 
exit condition is M 2 = 1, but the exit pressure is unknown. The details of the procedure 
are as follows: 

Choked film thickness evaluation . - When the flow is just choked, it is known that 
Pg = Pg and Mg = 1. By combining equation (26), the relation for Pq/Pj, with equa- 
tion (20), the equation for P^/Pg = (Pj/Pg) in the seal leakage passage, one obtains 



Since Pq/Pg is known, the Mach number M^ can then be determined from this equa- 
tion. Once the Mach number is known, all other flow quantities can be determined. The 
minimum choked film thickness is calculated as follows: 

(1) A film thickness h is assumed (chosen to be 1 mil in the computer program). 

(2) A Reynolds number is calculated from this film thickness. 

(3) This Reynolds number determines the flow regime, laminar or turbulent, and 
therefore the friction factor, (f = constant/Re 11 , where the constant and exponent n are 
program inputs and depend on whether the flow is laminar or turbulent. ) 

(4) A film thickness is then calculated by use of equation (16): 


B(M.) = — 

1 D 


2f AR 
h 


(29) 
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(For radial flow between coaxial parallel disks and parallel plates, the hydraulic diam- 
eter D = 2h. ) If this h does not agree with the previous h, steps (2) to (4) are repeated 
until the h's agree. This defines the critical or minimum choked film thickness h* . 

It should be noted that equation (28) cannot be solved for all values of Pg/Pg, C^, and 
y. Consequently, h* cannot be found. This condition arises for very low subsonic 
flows; hence, the flow is assumed to be subsonic. 

Sub sonic flow case . - For input film thickness less than the minimum choked film 
thickness, the flow is subsonic and the exit condition is Pg = Pg. In this case it is con- 
venient to use the length defined by equation (14); that is, refer all quantities to a 
fictitious length of the seal for which the condition at the exit is sonic (choked). To start 
the iteration procedure, the entrance Mach number is assumed to be the same as 
for critical flow. All other quantities at the inlet are then calculated. The value of the 
friction parameter B(M^) is found from equation (11) and the choking length from equa- 
tion (16); that is, Lj = B(Mj)D/4f. The length Lg can be calculated from Lg = Lj - AR 
and then the friction parameter B(Mg) from equation (11). Once B(Mg) is known, the 
Mach number Mg can be determined and, hence, the pressure Pg. This procedure is 
then repeated with a new Mach number Mj until Pg(M^) = Pg. 

Cho ked flow case . - For choked flow, the input film thickness is greater than the 
choking film thickness and the exit condition is Mg = 1. As in the subsonic flow case, 
the iteration procedure begins with the entrance Mach number assumed to be the choking 
film thickness entrance Mach number. Again all entrance quantities are calculated. The 
friction factor f is calculated from the Reynolds number, and the friction parameter 
B(Mj) = 4f L/D is calculated from equation (11). The trial flow length L is calculated 
as 


B(M-i )D 

L = (30) 

4f 

If L is not equal to the true flow length AR, the procedure is repeated with a new Mj 
until L does agree with AR. 

Calculation of Variables at Intermediate Points 

Once the entrance conditions are known, the Mach number, temperature, velocity, 
and pressure are determined at each desired location X N along the seal passage length 
in the following way: 

(1) L n is found from equation (14). 


13 



(2) B(M^) is calculated from equation (16). 

(3) is found by an iterative solution of equation (11). 

(4) T n is found from equation (19) where T 2 = T* and M 2 = 1; hence, 



(5) The velocity u^ is found from the Mach number definition. 

(6) The density is found from the mass conservation equation (1): 


P N = 


Pl u l 

U N 


(7) The pressure is found from the perfect gas law, equation (6). 


( 31 ) 


Computer Program 

The computer program in reference 12 used to carry out this analysis is written in 
FORTRAN IV. The program is described in detail including all input and program vari- 
ables and output options. Program listing and flow charts are presented and a sample 
problem with computer printout is given. The input and output can be in either the U.S. 
Customary or International System of Units. 


Remarks on Physics of Flow 

Static pressure drops to overcome flow friction. This pressure drop increases the 
specific volume of the fluid. Since the area change is negligible, the mean velocity 
must increase as the specific volume increases in order to maintain the same mass flow 
rate at each section of the leakage path. As the velocity increases the fluid momentum 
also increases, which requires an additional pressure drop. This results in a still 
greater increase in specific volume. This process will continue until the end of the 
leakage path or until the fluid reaches the maximum (choking) condition which will occur 
when the Mach number is one at the exit (exit velocity is sonic). Should the Mach num- 
ber reach a value of one somewhere along the leakage path interior, it is expected that 
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behavior similar to duct flow will occur. For duct flow it can be shown that the flow 
process will adjust itself until the point at which the Mach number reaches one is shifted 
to the exit of the leakage passage. The mass flow rate of the fluid is the maximum 
which can be handled for a given inlet density and passage cross-sectional area. This 
flow process is known classically as Fanno line flow. 

The entrance velocity head loss is negligible in subsonic viscous flow compared to 
the total seal pressure drop; however, under choked conditions the entrance velocity 
head loss is no longer negligible. This results in an entrance pressure loss. Also 
under choked flow conditions, the exit pressure is larger than the sump pressure and 
increases with film thickness. The pressure decreases through expansion waves to the 
sump pressure in the outer cavity. However, these expansion waves do not significantly 
affect the axial force balance on a seal nosepiece. 

As previously stated, in most face seals the area expansion is negligible due to the 
radius ratios being close to one. Hence, if the flow becomes choked it will usually be 
at the exit. Choking can occur at the entrance under some conditions. If there is flow 
separation at the entrance (due, e. g. , to the flow turning into the sealing faces), there 
could be a large entrance area decrease with choking occurring at the vena contracta. 
For smaller radius ratios such as those that characterize externally pressurized gas 
bearings, choking will occur at the entrance. 



RESULTS AND DISCUSSION 

The computer program from reference 12 was used to carry out the quasi- one- 
dimensional flow analysis. Several problems have been solved by using this program. 
One such problem is a case used in the design study of reference 5. The sealed pres- 
sure and temperature were 45 newtons per square centimeter absolute (65 psia) and 
311 K (100° F), respectively. The exit ambient pressure was 10.3 newtons per square 
centimeter absolute (15 psia). The sealing dam radial width (seal leakage passage 
length) was 0. 127 centimeter (50 mils). Figure 5 shows the calculated seal gas leakage 
as a function of film thickness, which ranges from 0. 003 to 0.0051 centimeter (0. 12 to 
2 mils). 

The analysis agrees with the differential analysis of reference 11, which is the 
classical compressible viscous analysis for subsonic viscous flow when the Mach number 
is less than 1/y/y. As shown in figure 5, for film thicknesses less than 8 micrometers 
(0. 3 mil), the Mach number is less than 1 /\/y = 0. 845 and the leakage dependence on 
film thickness is the classical cubic dependence. The isothermal viscous flow model 
loses its validity when the Mach number exceeds l/^y. The present approximate 
model, however, is valid for both subsonic and choked flow. Figure 5 shows that for 
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Figure 5. - Seal gas leakage as function of film thickness for parallel sealing dam surfaces. Radial dam width, AR, 0. 130 centimeter 
(50 mils); sealed air pressure, Pq, 45 newtons per square centimeter absolute (65 psia); sealed air temperature, Tq, 310 K (100° F); 
sump pressure, P 3 , 10 newtons per square centimeter absolute (15 psia). 


gaps larger than 8 micrometers (0. 3 mil) the leakage has a less than cubic dependence 
on film thickness. Note that choking occurs at a film thickness of 13 micrometers 
(0. 52 mils). The limiting case of orifice flow, in which the flow rate varies linearly 
with film thickness, would be achieved when the film thickness approaches the order of 
the sealing dam width of 0. 127 centimeter (50 mils). 

As indicated in figure 5, transition from laminar to turbulent flow occurs for a film 
thickness of approximately 0. 0033 centimeter (1. 3 mils). A Reynolds number, based on 
hydraulic diam eter, equal to 2300 has been chosen to be the critical transition Reynolds 
number. This appears to be a universal critical transition Reynolds number for flows 
in ducts, pipes, and bearings. For laminar flow, a mean Fanning friction factor of 
24/Re 2h was used. This friction factor is derived from the classical, viscous compres- 
sible flow solution in reference 11 and is shown in appendix D. For turbulent flow, the 
Blasius friction factor - Reynolds number relation was chosen. The mean friction 
factor equals 0. 079/Re*^. For the transition flow regime the exact nature of the flow 
(hence, friction factor) is complex and not fully understood. The friction factor used 
for this flow regime is derived in reference 12 by assuming a smooth transition from 
the laminar friction factor to the turbulent friction factor. Reynolds numbers in the 
2300 to 3000 range are arbitrarily selected as the transition flow regime. 
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Correlation With Experiment 


Experiments were conducted at the Lewis Research Center to check the validity of 
the model for choked flow and for a radial geometry representative of face seals. 

Leakage flow was studied for a face seal geometry (14 cm (5. 50 in. ) inner diameter 
15 . 2 cm (6. 00 in. ) outer diameter, separated by a fixed parallel gap of 0. 0038 cm 
(1. 5 mils)). The reservoir pressure of 41. 8 newtons per square centimeter absolute 
(60. 6 psia) was held fixed while the exit pressure varied. Figure 6 shows a comparison 
of this analysis, the classical viscous subsonic flow differential analysis, and the radial 
flow experiment. Notice that the experiment shows that the flow does become choked. 
The maximum flow rate (choked) is 0. 0080 kilogram per second (0. 0176 lbm/sec). The 
present analysis predicts a slightly higher mass flow but agreement is within 19 percent 
over most of the range. On the other hand, classical compressible viscous flow theory 
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Figure 6 . - Comparison of classical viscous flow theory and present analysis with 
radial Now experimental results. Sealed air pressure, Pq, 41. 8 newtons per 
square centimeter (60.6 psia); sealed air temperature, Tq, 300 K ( 8 C P F ); film 
thickness, h, 0.0038 centimeter (1. 5 mils); entrance radial width, Rj, 6.98 
centimeters (2. 75 in. ); exit radial width, R 2 , 7.62 centimeters (3.00 in. ); sump 
pressure P 3 varied. 
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overestimates the flow rate considerably at all pressure ratios except very near one. 

The classical theory predicts no limiting mass flow. However, if the limit is chosen as 
the -mass flow when the Mach number is unity, predicted flows are about 80 percent 
higher than those observed experimentally for choked flow. 

The present analysis used a friction factor of 24/1162^ over the laminar flow region. 
The good agreement in figure 6 indicates that this value should be sufficient for engineer- 
ing purposes. The important point is that it was derived analytically from classical 
compressible viscous flow theory. As is the case with duct flow, it appears to be nearly 
invariant with compressibility. Exact agreement could have been achieved if a different 

friction factor had been used. In this case almost exact agreement with the experimental 

— n 

curve would be achieved if f were chosen to be 6/Re ; however, it may not be a 
reliable predictor for other experiments and applications. 

There are many possible reasons for the disagreement between the analytical and 
experimental results. It is extremely difficult to maintain a uniform film thickness for 
the gap sizes of interest as small as 0. 0038 centimeter (1. 5 mils). Problems in clamp- 
ing, machining tolerances, and distortions due to large pressure differentials result in 
experimental error. Another possible source of disagreement between theory and exper- 
iment could be due to area expansion, which the theory neglects. In the experiment, the 
inner diameter to outer diameter ratio is 0. 92. Thus one could expect a lower flow rate 
than calculated. (Calculated mass flow is based on the arithmetic mean diameter. ) 

The theory also assumed isentropic entrance flow conditions. For large pressure 
differentials the entrance velocity is large enough to cause a substantial entrance pres- 
sure drop. The result of accounting for this entrance loss condition would also result 
in a decrease in mass flow. As illustrated in figure 6, an entrance velocity loss coeffi- 
cient equal to 0. 6 gave excellent agreement with experiment. 

Since the Mach number is defined using the mean flow velocity, the flow will choke 
locally at the centerline sooner than that predicted here at unit Mach number. The cri- 
teria used in this analysis would be accurate if a slug profile were present. Therefore, 
the prediction may be better for turbulent flow. Because physically there is no slip on 
the walls, the flow should never have an exact limit. The experimental data in figure 6 
show this; however, as shown in figure 6, as the pressure ratio is decreased 

below 0.25, the mass flow can be considered choked for engineering purposes. Note the 
asymptotic behavior which shows a very small increase in mass flow with decreasing exit 
pressure. 

The radial flow experiment was repeated for two other reservoir pressures. The 
experimental results are shown in figures 7 and 8 and compared with the analysis. For 
a reservoir (sealed) pressure of 27.6 newtons per square centimeter absolute (40 psia), 
the flow was laminar over the entire range studied (fig. 7). Agreement with the analy- 
sis, assuming isentropic entrance conditions, was within 18 percent. With an entrance 
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Figure 7. - Comparison of present analysis with radial flow experimental 
results. Sealed air pressure, Pq, 27.6 newtons per square centimeter 
absolute (40. 0 psia); sealed air temperature, To, 296 K (73° F) ; film 
thickness, h, 0.0038 centimeter (1.5 mils); entrance radial width, R 1( 
6.98 centimeters (2.75 in.); exit radial width, R 2 , 7.62 centimeters 
( 3. 00 in. ) ; laminar flow exists everywhere. 
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Figure 8 . - Comparison of present analysis using Blasius friction factor (f = 0.079/ 
Re 1/4 ) with radial flow experimental results. Sealed air pressure, Pq, 67. 1 new- 
tons per square centimeter absolute (97. 3 psia); sealed air temperature, Tq, 

295 K (71° F); gap, 0. 0039 centimeter (1.53 mils); entrance radial width, R 1( 

6. 98 centimeters (2. 75 in. ); exit radial width, Ify 7. 62 centimeters (3. 00 in. ). 



loss coefficient of 0. 6, agreement between analysis and experiment was within 5 percent. 

Figure 8 shows the 67. 1 newtons per square centimeter absolute (97.3 psia) reser- 
voir pressure case experimental results compared with the analysis using the Blasius 
turbulent mean friction factor of 0.079/Re ' . The inverse one- fourth dependence on 

Reynolds number by the friction factor occurs for many fully developed turbulent flows. 
The analysis with isentropic entrance conditions overestimates the flow by 13 percent; 
however, with an entrance loss coefficient of 0. 6, there is excellent agreement with the 
experimental data. In the laminar flow regime the agreement is also excellent. Note 
the transition Reynolds number of 2300 occurred for a pressure ratio of 0. 87. The tur- 
bulent flow regime was set to begin at a Reynolds number of 3000 which occurred at a 
pressure ratio of 0. 55. 

Since shaft face seals are pressure balanced, it is necessary to know the sealing 
dam opening force which is found from 


F = W 


/ 


R 0 -R, 


(P- 


P 3 ) dX 


(32) 


For classical viscous, compressible flow, equation (32) is easily integrated in closed 
form to yield 


F = 


2W(R 2 - R^Pj 


1 - 




- P 3 (R 2 - R X )W 


(33) 


Note that the opening force is independent of the fluid properties and film thickness. 

This gives the well known result that parallel surfaces yield no film stiffness. However, 
this is not the case when the classical compressible viscous flow theory is no longer 
valid. Figure 9 shows a plot of seal opening force as ,i function of film thickness ob- 
tained using the computer program in reference 12 for both isentropic and 0. 6 entrance 
loss cases. Notice that for small film thickness (corresponding to M < 1 /\fy), the 
force is constant, as predicted by classical analysis. However, as film thickness (and 
Mach number) increases, the force actually increases slightly for the isentropic en- 
trance case. This is a condition of negative film stiffness. It would be undesirable to 
operate in this region unless an auxiliary film stiffness generating device, such as self- 
acting lift pads, could make the overall seal stiffness positive. This negative stiffness 
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Figure 9. - Sealing dam opening force for parallel surfaces comparing classical theory with 
present analysis. Sealed air pressure, Pq, 148 newtons per square centimeter absolute 
(215 psia); sump pressure, P 3 , 10 . 3 newtons per square centimeter absolute (15 psia); 
sealed air temperature, Tn, 700 K (800° F); sealing dam radial width, AR, 0. 127 centimeter 
(50 mils). 


may be one of the reasons why conventional gas film face seals T ’chatter Tf under certain 
high pressure conditions. 

Also notice in figure 9 that the opening force attains a peak value and then sharply 
decreases with increasing film thickness. This sharp decrease indicates a high positive 
film stiffness which heretofore was not thought to be present for parallel surfaces. 

The results using the analysis with an entrance loss coefficient of 0. 6 are also 
shown in figure 9. The negative stiffness region has been greatly reduced; however, 
there is a larger region of positive film stiffness. A maximum film stiffness of about 
102 000 newtons per centimeter (57 800 lbf/in. ) occurs at a film thickness of 0. 00152 
centimeter (0. 6 mil). Notice that there is another possible negative stiffness region 
where the flow is in transition between laminar and turbulent flow. The uncertainty of 
the real flow behavior in this region may mean different behavior than that calculated. 
Positive film stiffness also occurs in the turbulent flow regime. Of course, in an actual 
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seal design, the advantage gained by this added film stiffness must be weighed against 
the higher leakage flow resulting from choked flow operation. 


Limitations of the Analysis 

One of the limitations of this analysis is the neglect of the details of the flow in the 
entrance length. It is common in lubrication theory to neglect the entrance region en- 
tirely. This is reasonable for the slow viscous flows which characterize lubrication 
flows. However, when the flow becomes choked, the entrance velocities are high and it 
is possible that this neglect is improper. 

The model developed herein may still be applied with good accuracy in cases where 
entrance effects are significant by using an entrance loss coefficient (see eqs. (25) to 
(27) and appendix C). Fleming and Sparrow have shown in reference 16 the bulk of the 
entrance losses occur in a small region very close to the duct inlet. Laminar incom- 
pressible entrance loss coefficients have been calculated and measured for a variety of 
duct shapes (e.g. , refs. 16 and 17). Turbulent loss coefficients are not so widely re- 
ported; however, they may be calculated for parallel plate channels from the information 
in reference 18. The values for head loss coefficients k* are generally less than 20 
percent of those for laminar flow. Incompressible loss coefficients may be adequate for 
use with compressible flow (see appendix C). 

The use of the quasi- one- dimensional model can be justified by the good agreement 
with experiments. Also, calculations can be made very rapidly on a digital computer. 


Seal Design Example 

NASA has designed a seal with self-acting lift pads for potential use as a mainshaft 
seal in advanced gas turbine engines. This seal is described in references 4, 5, and 
19. The sealing dam portion of this mainshaft face seal was designed and studied using 
the analysis presented herein. 

Figure 10 shows the effect of sealing dam radial width on leakage flow film thick- 
nesses from 2. 5 to 0. 25 micrometer (0. 1 to 1 mil). The four radial dam widths used in 
the study were 0.051, 0. 127, and 0.254 centimeter (5, 20, 50, and 100 mils). These 
plots illustrate that, from strictly a leakage point of view, the longest leakage path pos- 
sible is most desirable. However, it is shown in reference 5 that the leakage path length 
must be compromised from a force balance point of view when surface deformations 
occur. Also shown in figure 10 is the leakage flow rate - film thickness relation for a 
0o 001 centimeter (0. 4 mil) sealing dam radial width. Notice that at film thicknesses 
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Figure 10. - Effect of seal face radial width on seal leakage. Sealed air pressure, Pg, 148 newtons per square centimeter absolute 
(215 psia); sump pressure, P 3 , 10. 3 newtons per square centimeter absolute (15 psia); sealed gas temperature, 700 K <800P F); 
parallel sealing faces. 
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greater than 0.00063 centimeter (0.25 mil) the leakage coincides with the leakage values 
obtained for a knife edge using the theoretical orifice flow equation. Also notice the 
linear variation of leakage flow with film thickness increase. For this narrow sealing 
dam width case when flow behaves as knife edge flow, the velocity entrance loss coeffi- 
cient is the same as the flow discharge coefficient used in orifice flow analysis. 

Figure 11 shows mass leakage flow as a function of pressure ratio for a fixed gap of 
0. 0001 centimeter (0 Q 4 mil) and varying sealing dam width. Notice the knife edge, as 
expected, has the largest leakage flow rate. These values have been calculated from 
the theoretical orifice flow equation. The computer program result for a sealing dam 
width of 0. 0001 centimeter (0. 4 mil) agreed with the theoretical orifice flow equation. 

As the sealing dam width increases, the critical pressure ratio for choking decreases 
as indicated by the sonic line in figure 11. Again the advantage of wider (longer) sealing 
dams is apparent. 

Recent NASA sponsored tests (ref. 19) of a shaft face seal with self-acting lift aug- 
mentation have demonstrated the feasibility of operation at gas temperatures up to 910 K 
(1200° F), pressure differentials across the seal up to 172 newtons per square centime- 
ter (250 psi), and relative surface speeds up to 130 meters per second (450 ft/sec). 
These tests simulated the engine operating conditions such as takeoff, climb, cruise,, 
and descent. The seal appears to be operating as predicted by the design analysis. 
Figure 12 shows the nosepiece seal assembly after the 120- hour endurance test. Fig- 
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Figure 13. - Closeup of carbon nosepiece of gas-film seal after 200 hours of endurance testing. Total 
time on seal, 338.5 hours. 


ure 13 shows a closeup of the carbon nosepiece from the shaft face seal with self-acting 
lift pads after 320 hours of steady- state endurance testing. The total time of testing on 
the carbon nosepiece was 338. 5 hours. The pretesting lapping marks are still observ- 
able. A surface profile trace indicated that the average wear on the surface was less 
than 0. 13 micrometer (5 p. in . ). A few shallow scratches (1 to 2 jum, 50 to 100 jUin. ) 
were noticed. The test was concluded after 500 total hours of carbon nosepiece opera- 
tion. The seal faces had encountered over 50 startups and shutdowns. The leakage 

_ O O 

rates varied from 5. 19x10 to 15. 6x10 standard cubic meter per second (11 to 32 std 

ft /min) with leakage generally averaging about 1. 18x10 ° standard cubic meter per 
3 

second (25 std ft* /min). These leakage rates were close to those theoretically predicted. 


SUMMARY OF RESULTS 

An analysis has been presented for compressible fluid flow across shaft face seals. 
This quasi- one-dimensional integral analysis includes fluid inertia and entrance losses 
in addition to viscous friction which is accounted for by a mean friction factor. Subsonic 
and choked flow conditions can be predicted and analyzed. The model is valid for both 
laminar and turbulent flows. The following pertinent results were found: 

1. Results agree with the classical subsonic compressible viscous flow theory for 
Mach numbers less than l/y^y. Excellent agreement with experiment is achieved if an 
entrance velocity loss coefficient of 0. 6 is used in the laminar flow regime. The Blasius 
friction factor of 0. 079/Re in the analysis yielded good experimental agreement in 
the turbulent flow regime. 
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2. Near and at critical flow conditions (Mach number = 1) the isentropic entrance 
flow analysis shows a negative film stiffness; whereas for choked flow the analysis pre- 
dicts: a high positive film stiffness. The analysis with an entrance loss coefficient of 
0. 6 predicts a high positive film stiffness for choked flow and yielded a maximum posi- 
tive film stiffness of 102 000 newtons per centimeter (57 800 lbf/in. ) for a design prob- 
lem-studied. The choked flow results contrast to classical compressible viscous flow 
results which show no film stiffness for parallel surfaces. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 19, 1972, 

132-15. 
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APPENDIX A 


SYMBOLS 


A gap cross-sectional area 

B friction parameter, 4f L/D 

C L velocity entrance loss coefficient 

Cp specific heat at constant pressure 
D hydraulic diameter, 2h 

e specific internal energy 

F sealing dam force 

f Fanning friction factor, t /(pu /2) 

f mean Fanning friction factor 

h film thickness (gap) 

i specific enthalpy 

k' head loss coefficient 

flow length from entrance to point of choking 
M Mach number, u/^ydlT 

M mass flow 

P pressure 

Q volume leakage flow rate 

r radius 

AR sealing dam radial width (physical flow length), Rg - R-^ 
<R gas constant = universal gas constant/molecular weight 
Re Reynolds number, pu2h/p 

T temperature 

u average velocity 

W flow width 

X radial coordinate direction 

z coordinate across film thickness 



p absolute (dynamic) viscosity 
y specific heat ratio 
p density 
r shear stress 
Subscripts: 

N location along flow leakage length 
w wetted surface 

0 sealed (reservoir) conditions 

1 entrance conditions 

2 exit conditions 

3 ambient sump conditions 
Superscript: 

~ average value 
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APPENDIX B 


FRICTION FACTOR - MACH NUMBER RELATION 


Equation (9) may be derived in the following way. From the equation of state 
(eq. (7)) and the energy equation (eq. (5)) one obtains 


dP = dp + dT _ dp /y-y du 2 
P p T p \ 2 / u 2 


(Bl) 


Combining this with the mass conservation equation (eq. (2)) gives 


dP 1 + (y - 1)M 2 
P 2 



(B2) 


Combining equation (B2) with the momentum equation (eq. (8)) yields 


2yM 2 f 

D 



(B3) 


2 2 

From the definition of the Mach number (M = u /y(RT) and the energy equation, one 
obtains 


V 2 / u 2 



Equations (B3) and (B4) then give 


2 yM 2f dx 

D 


(1 - M 2 ) dM 2 

2i + (y 1 ) M 2 M 2 
2 


(B4) 


(B5) 
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from which equation (9) can be obtained; that is, 


APPENDIX C 


REMARK ON LOSS COEFFICIENT 


There are several effects which cause nonisentropic pressure drops in the face seal 
entrance region. Because of an abrupt geometric change at the seal passage entrance, 
there is a nonuniform profile caused by flow turning, separated flow, and ’’vena con- 
tracta’’ effect which results in an entrance pressure drop. The pressure drop in the 
flow- development length is higher than that in the fully developed flow region because of 
two effects (ref. 16). The first effect is higher wall shear caused by higher transverse 
velocity gradients. The second effect is the momentum increase as the velocity distri- 
bution becomes less uniform. 

The role of the entrance velocity loss coefficient may be better understood by con- 
sidering the incompressible Bernoulli equation relating the stagnation reservoir pres- 
sure with the static and dynamic pressure at the seal entrance: 

P 0 = p i + ipi(ui) 2 (Cl) 

u 1 2 1 'ideal 


The actual velocity at the seal entrance is less than ideal; it can be expressed as 


Ui = C n 


H 


ideal 


(C2) 


where C L is the entrance velocity loss coefficient (see ref. 20). Substitution of equa- 
tion (C2) in (Cl) and use of the Mach number definition yield 


P 0 - P 1 + 


yPiM 2 

2C? 


or 


p i = 


i + 


yM 2 

2C 2 


(03) 


(C4) 
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A binomial expansion of the denominator of the compressible entrance pressure 
equation (26) yields 


1 1 + - (r - l) — 


]y/(r-l) 


= 1 +- 

2 


1 yM 2 | 1 yM 4 


8 r*,4 


y( 2 - y) M 
48 


6 


(C5) 


The first two terms on the right side strongly predominate for Mach numbers less 
than unity. For example, when M/C^ = 1 and y = 1. 4, 


i 1 / oM 2 
1 + - (y - 1) 

2 r 2 

h 

The error is less than 11 percent if only the first two terms are used. The accuracy 
would be greater at lower Mach numbers. The first two terms of the expansion are the 
same as the denominator of the Bernoulli equation (C4) for a gas behaving as a quasi- 
incompressible fluid. Using this observation the relation between the entrance velocity 
loss coefficient C-^ and head loss coefficient k' commonly used in hydraulics can be 
found. The incompressible Bernoulli equation with a head loss coefficient is 


y/(y-l) 


= 1 + 0.70 + 0. 18 + 0. 02 + 


p U 2 


p l u l 


+ Pi 


or 


P 


1 



(k* + 1) 


Pi u j 

2 



- (k* + 1) 


y p i M j 

2 


Comparing equation (C7) with equation (C3) yields 



(C6) 


(C7) 


(C8) 
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This relation enables one to convert head loss coefficients reported in the literature 
to velocity loss coefficients C^, which are more convenient to use in Mach number 
relations. 


APPENDIX D 


DERIVATION AND JUSTIFICATION 
Derivation of Mean Fanning Friction Factor 

The mean Fanning friction factor commonly used in fluid mechanics is defined as 

I = D dP 

2pu 2dx 

The classical viscous radial velocity distribution is (ref. 11) 

u = -L ^ ( Z 2 _ hz) 

2/i dx 


The mean radial velocity is found from 



h 2 dP 
12 p dx 


or 

_ dP _ 12 pu 
dx h 2 

Since D = 2h, the mean Fanning friction factor becomes 

f =21 
Re 

Note the Fanning friction factor is one- fourth the Darcy- Weisbach friction factor which 
is commonly used in hydraulics. (In the text the velocity u is understood to be the 
mean velocity. ) 
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Validity of Mean Friction Factor Evaluation Using Entrance Reynolds Number 

The question arises if it is proper to use the entrance Reynolds number to evaluate 
the mean friction factor. Examining the Reynolds number, 

Re =^H2h 

M 


or 


Re - 2M _ constant 
Wju M 

Hence Re varies inversely with the absolute viscosity of the fluid which in turn varies 
with the static temperature. 

To see the maximum possible change, the maximum temperature possible occurs 
when = 0 and M2 = 1. 

Thus, the exit- entrance temperature ratio can be found using equation (19) and is 



1 + 


1 + 


y - 1 
2 

y - 1 
2 



2 

y + 1 


For y = 1.4, Tg/T. = 5/6. This represents about a 17 percent variation in the static 
temperature for the worst case. From Sutherland’s law 


11 cc T 


0.76 


f cc Re - cc jU oc T®' 

Hence, the maximum viscosity change along the seal leakage path will be small, and the 
entrance conditions can be used satisfactorily to evaluate the mean friction factor. 
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